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We reelaborate on the basic properties of lossless multilayers. We show that the transfer matrices 
for these multilayers have essentially the same algebraic properties as the Lorentz group S0(2,l) in 
a (2 + l)-dimensional spacetime, as well as the group SL(2,R) underlying the structure of the ABCD 
law in geometrical optics. By resorting to the Iwasawa decomposition, we represent the action of 
any multilayer as the product of three matrices of simple interpretation. This group-theoretical 
structure allows us to introduce bilinear transformations in the complex plane. The concept of 
multilayer transfer function naturally emerges and its corresponding properties in the unit disc are 
studied. We show that the Iwasawa decomposition reflects at this geometrical level in three simple 
actions that can be considered the basic pieces for a deeper undestanding of the multilayer behavior. 
We use the method to analyze in detail a simple practical example. 

PACS numbers: 120.5700 Reflection, 120.7000 Transmission, 230.4170 Multilayers, 000.3860 Mathematical 
methods in physics 



I. INTRODUCTION 

Layered media are important for many applications in 
modern optics. To fully capitalize these media one must 
have a clear picture of all the mechanisms involved in 
the propagation of the optical waves in them. In spite of 
this fact, multilayer optics is usually approached from a 
practical view, in relation to optical filters and the like. 
In this spirit, the topics covered in most of the textbooks 
on the subject use a mixture of design, manufacture, and 
applications; treating only the basic physics needed to 
carry out practical computations |l[]. 

However, for a variety of reasons layered media have 
physical relevance on their own ||, ^ . The mathematical 
basis for understanding their behavior is the fact that 
the matrix representing any lossless multilayer belongs 
to the group SU(1,1). It is known that this group is lo- 
cally isomorphic to the (2-1- l)-dimensional Lorentz group 
S0(2,l) 1^, This leads to a natural and complete iden- 
tification between reflection and transmission coefficients 
of the multilayer and the parameters of the correspond- 
ing Lorentz transformation. It is, precisely, the abstract 
composition law of SU(1,1) the ultimate responsible for 
the curious composition law of these reflection and trans- 
mission coefficients. Moreover, this fact allows one to 
perform experimental tests of special relativity with sim- 
ple optical measurements involving multilayers 0, 

In this respect, another remarkable fact to be consid- 
ered is that SU(1,1) is also isomorphic to SL(2,IR), which 
is the natural group underlying the mathematical struc- 
ture of the celebrated ABCD law in first-order optics. In 
this paper, we also exploit this correspondence to explore 



intriguing connections of layered media with geometrical 
optics. 

These purely algebraic results seem to call for a geo- 
metrical interpretation. It is difficult to overestimate the 
role played by geometrical ideas in all the branches of 
physics, particularly in special relativity. The hyperbolic 
geometry associated with the group SU(1,1) [or, equiva- 
lently, S0(2,l)] is an estabfished doctrine [||. In fight of 
these considerations, it is easy to convince oneself that 
this geometrical approach might provide deeper insights 
into the action of a multilayer in a wider unifying frame- 
work that can put forward fruitful analogies with other 
physical phenomena. 

In consequence, it is natural to view the action of a 
SU(1,1) multilayer matrix as a bilinear (or Mobius) trans- 
formation on the unit disc, obtained by stereographic 
projection of the unit hyperboloid of S0(2,l). This kind 
of bilinear representations have been discussed in detail 
for the Poincare sphere in polarization optics , for 
Gaussian beam propagation |p^ , and are also useful in 
laser mode- locking and optical pulse transmission p^ . 

In addition, the isomorphism with SL(2,R) allows us 
to translate the geometrical structure defined in the unit 
disc to the complex plane, recovering in this way an al- 
ternative model of the hyperbolic geometry that is useful 
in shedding light into the behavior of the multilayer. 

In spite of these achievements, the geometrical action 
of an arbitrary lossless multilayer could still become cum- 
bersome to interpret in physical terms. In fact, in prac- 
tice it is usual to work directly with the numerical values 
of a matrix obtained from the experiment, which cannot 
be directly related to the inner multilayer structure. To 
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remedy this situation, we have resorted recently to 
the Iwasawa decomposition, which provides a remarkable 
factorization of the matrix representing any multilayer 
(no matter how complicated it could be) as the product 
of three matrices of simple interpretation. 

At the geometric level, such a decomposition trans- 
lates directly into the classification of three basic actions, 
which are studied in this paper, that are the basic bricks 
from which any multilayer action is built. 

The contents of this paper are organized as follows. 
In Section 2 we present in detail some relevant algebraic 
properties of lossless multilayers. These properties have 
a direct translation into hyperbolic geometry, both in the 
unit disc and in the complex plane, which is explored in 
Section 3. Finally, Section 4 devote to develop a simple 
yet relevant example that illustrates the power of this 
approach and to expose our concluding remarks. 



II. SOME PROPERTIES OF LOSSLESS 
MULTILAYER MATRICES 

A. Transfer matrix for a lossless multilayer 

We first briefly summarize the essential ingredients of 
multilayer optics we shall need for our purposes [ p^ . 
The configuration is a stratified structure, illustrated in 
Fig. 1, that consists of a stack of 1, . . . , j, . . . , m, plane- 
parallel lossless layers sandwiched between two semi- 
infinite ambient (a) and substrate (s) media, which we 
shall assume to be identical, since this is the common 
experimental case. Hereafter all the media are supposed 
to be lossless, linear, homogeneous, and isotropic. 

We consider an incident monochromatic linearly polar- 
ized plane wave from the ambient, which makes an angle 
9q with the normal to the first interface and has ampli- 
tude Ei'^K The electric field is either in the plane of 
incidence {p polarization) or perpendicular to the plane 
of incidence (s polarization). We consider as well another 
plane wave of the same frequency and polarization, and 

with amplitude Es ■* , incident from the substrate at the 
same angle 9o- 

As a result of multiple reflections in all the interfaces, 
we have a backward-traveling plane wave in the ambi- 
ent, denoted \ and a forward-traveling plane wave 



in the substrate, denoted E', 



(+) 



If we consider the field 



amplitudes as a vector of the form 



E 



Ei-) 



(1) 



which applies to both ambient and substrate media, then 
the amplitudes at each side of the multilayer are related 
by a 2 X 2 complex matrix Mqs, we shall call the multilayer 
transfer matrix pq], in the form 



Mas = 





a 






(3* 


a 



(3) 



where the complex numbers 
Ras = \Ras\ exp(ip). 



Tas = \Tas\eMiT), (4) 



are, respectively, the overall reflection and transmission 
coefficients for a wave incident from the ambient. Note 
that 



1 - \Ras 
\Tas? 



= +1. 



(5) 



Therefore, the condition det Mas ~ +1 is equivalent to 
l^asP + iT'asP = Ij and then the set of lossless multilayer 
matrices reduces to the group SU(1,1), whose elements 
depend on three independent real parameters. 

When the ambient and substrate media are different, 
this result also holds after a convenient renormalization of 
the field amplitudes Q . The identity matrix corresponds 
to Tas = 1 and Ras — 0, so it represents an antireflec- 
tion system. The matrix that describe the overall system 
obtained by putting two multilayers together is the prod- 
uct of the matrices representing each one of them, taken 
in the appropriate order. So, two multilayers, which are 
inverse, when composed give an antirefiection system. 



B. A basic factorization for multilayers: the 
Iwasawa decomposition 

Many types of matrix factorizations have been consid- 
ered in the literature |9j , all of them decomposing 
the matrix as a unique product of other matrices of sim- 
pler interpretation. Particularly, given the essential role 
played by the Iwasawa decomposition, both in fundamen- 
tal studies and in applications to several fields (especially 
in optics), one is tempted to investigate also its role in 
multilayer optics. 

Without embarking us in mathematical subtleties, the 
Iwasawa decomposition is established as follows ^ : 
any element 5 of a (noncompact semi-simple) Lie group 
can be written as an ordered product of three elements, 
taken one each from a maximal compact subgroup K, 
a maximal Abelian subgroup A, and a maximal nilpo- 
tent subgroup N. Furthermore, such a decomposition is 
global (in the sense that it applies to every group ele- 
ment) and essentially unique (in the sense that the el- 
ements of the factorization are uniquely determined in 
terms oi g). 

For the problem at hand of a lossless multilayer matrix 



e SU(1,1), the decomposition reads as |14 
Ma,s = K(0)A(e)N(i.) , 



(6) 



where 



(2) 



exp(i0/2) 






exp(— 10/2) 
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A(0 = 
N(z.) - 



cosh(^/2) i sinh(^/2) 
-i sinh(^/2) cosh(^/2) 

l~iiy/2 vjl 
v/2 l + iv/2 



(7) 



The parameters 0, ^, and v are given in terms of the 
elements of the muhilayer matrix by 



0/2 = arg(a + i(3) , 
i/2 ln(l/|a + j/3|) , 
1^/2 = Re{a(3*)/\a + iP\^ 



(8) 



where the ranges of the parameters are ^, G R and 
— 27r < (j) < 2tt. Therefore, given a priori hmits on a and 
fi (i.e., on Tas and i?as), one could easily establish the 
corresponding limits on the parameters 0, ^ and v, and 
vice versa. 

All the matrices in Eq. (||) are of SU(1,1) and, there- 
fore, leave invariant the expression — at 
each side of the multilayer; i.e., 

\Ei+^' -\Ei-^? ^\Ei-^^' ~\Ei-^\\ (9) 

which is nothing but the energy-flux conservation, as 
could be expected from physical considerations. In ad- 
dition, the matrix K{(j>) preserves the product E'^^^E^~\ 

the matrix A(^) preserves the quadratic form + 

, and the matrix N(i^) preserves the sum -|- 
iE(-l 



C. Connections with geometrical optics 

It this subsection we wish to show a remarkable anal- 
ogy with geometrical optics that could provide new in- 
sights into the interpretation of the action of any lossless 
multilayer. To this end wc note that by applying the 
unitary matrix 



U = 



1 
V2 



1 i 
i 1 



to both sides of Eq. (g) , we can recast it as 

where the new field vectors are defined as 

f£''-'\^U^^±(E^^\+^E^-^ 



e = 



and the conjugate multilayer matrix is 



a b 
c d 



(10) 

(11) 
(12) 

(13) 



where Mas is a matrix with dct Mas = +1 and whose 
elements are real numbers given by 



a = Rc(a) + Im(/3) , 
; = -Im(a) + Re{(3) , 



b = Im(a) + Re(/3) , 



Re(a) - Im(/3) 



(14) 



In other words, the matrices Mas belong to the group 
SL(2,]R), which plays an essential role in a variety of 
branches in optics, ranging from the celebrated ABCD 

to squeezed 
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law of geometrical optics [p2| , ^ ^ ^ ^ 

states in quantum optics ]26|, ^Tf, including beam prop- 
agation problems ||28|] . The transformation by U estab- 
lishes a one-to-one map between the group SL(2,M) of 
matrices Mas and the group SU(1,1) of matrices Mq^, 
which allows for a direct translation of the properties 
from one to the other. Such a correspondence is an iso- 
morphism of groups, because 



(15) 

By conjugating with U the Iwasawa decomposition m), 
we get the corresponding one for SL(2,M), which has been 
previously worked out p9[ : 



Mas=IC{cj))A{Om'^) , 



(16) 



where 



/C(0) = 
AiO = 



cos(0/2) sin(0/2) 
-sin(0/2) cos((?!)/2) 

exp(?/2) 
exp(-e/2) 

1 0" 
V 1 



(17) 



Now, we can interpret the physical action of the ma- 
trices appearing in both factorizations in SU(1,1) and 
SL(2,]R), respectively. In this way, K(0) represents the 
free propagation of the fields E in the ambient medium 
through an optical phase thickness of 0/2. Obviously, 
this reduces to a mere shift of the origin of phases. Al- 
ternatively, one can consider /C(0) as an ABCD matrix 
in geometrical optics that applies to position x and mo- 
mentum p (direction) coordinates of a ray in a transverse 
plane Ig^l- These are the natural phase-space variables 
of ray optics and then /C(0) would represent a rotation in 
these variables In the multilayer picture, f^"''' can 
be seen as the variable x and £'^~^ can be seen as the 
corresponding p. 

In Eq. (^, the second matrix A(^) represents a sym- 
metric multilayer (i.e., the reflection coefficient is the 
same whether light is incident on one side or on the op- 
posite side of the multilayer, and so ta — pa = i''i'/2) 
with transmission and reflection phase shifts of ta = 
and Pa = ±7''/2, and a transmission coefficient T/\ = 
sech(^/2). There are many ways to get this performance, 
perhaps the simplest one is a Fabry-Perot system com- 
posed by two identical plates separated by a transpar- 
ent spacer. By adjusting the refractive indices and the 
thicknesses of the media one can always get the desired 
values. Viewed in SL(2,E), A{i) represents a magnifier 
that scales x up to the factor m = exp(^/2) and p down 
by the same factor p^ . 

The third matrix, N(t/), represents a system having 
T|\| = cos(tn) exp(iri\i) and = sin(ri\i) exp(iri\i), with 
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tan(TN) = 1//2. The simplest way to accomplish this 
task is by an asymmetrical two-layer system. Using the 
analogy with the ABCD matrix in geometrical optics, 
A/'(j^) represents the action of a lens of power v |gg . 

Finally, in complete equivalence with the invariants 
found for the Iwasawa decomposition of SU(1,1), the ma- 
trix /C preserves the sum f^^^^ -I- f , A preserves the 
product and Af preserves £^'^\ In addition, 

the energy-flux conservation Eq. (|^) can be recast now 
as the invariance 



D. Connections with special relativity 

Apart from the geometrical-optics perspective devel- 
oped in the previous subsection, multilayer action can be 
viewed in a relativisticlike framework, that has proved to 
be very appropriate to understand some peculiarities of 
multilayer behavior. To this end, let us first recall some 
well-known facts about the Lorentz transformations in 
(2-|-l)-dimensions. 

Introducing a three-dimensional real vector space of 
vectors with components (x", x^, x^), where = ct, a 
Lorentz transformation A is a linear transformation be- 
tween two coordinate frames 



= A^ 



(18) 



(the Greek indices run from to 2), such that the pseudo- 
Euclidean bilinear form 



{x\y) = x"y° ~ xV - ^'y^ (19) 
^'\y') — (2^1^)7 which immedi- 



remains invariant; i.e 
ately implies 

detA^il. (20) 

Thus, the set of Lorentz transformations can be classified 
in two classes: proper transformations, with det A = +1, 
and improper ones, with det A = — 1. The proper Lorentz 
transformations form a subgroup, but the improper ones 
do not. 

Furthermore, it is easy to check that jA" o| > 1- There- 
fore, the transformations of the Lorentz group can also be 
classified according to the sign of A° the orthochronous 
Lorentz transformations, with Aq > 1, form a subgroup, 
but the antichronous ones, with A" < do not. 

We are interested in dealing with the three-parameter 
restricted Lorentz group S0(2,l); i.e., the group of the 
Lorentz transformations with determinant -1-1 and that 
do not reverse the direction of time. In fact, we wish 
to discuss now a very close correspondence between 
the group SU(1,1) introduced above and the restricted 
Lorentz group S0(2,l). We shall show this important 
correspondence ex plic itly in the following form that we 
recall for clarity [Q: with each point of our three- 
dimensional vector space with coordinates x^ we asso- 
ciate the Hermitian matrix 



X = 



I X 



I X 



(21) 



where ao = I is the identity, and cri and cr2 are the cor- 
responding Pauli matrices. Note, that detX = {x\x) = 

(x0)2-(xl)2-(x2)2. 

Now, if M e SU(1,1), then the matrix 



X' = M X M 



(22) 



where the symbol f denotes the Hermitian conjugate, in- 
duces a Lorentz transformation on the coefficients x^. 

It is clear from Eq. (^2|) that the matrices M and — M 
generate the same A, so this homomorphism is two-to- 
one. This equation can be easily solved to obtain M 
(uniquely defined up to the sign) from a given A. In 
fact, it is easy to find it explicitly as [Eol 



(23) 



While M acts on two-dimensional complex vectors like 
those in Eq. (1), the induced Lorentz transformation 
A(M) acts on three-dimensional real vectors of the form 
[which are the space-time counterparts of Eq. (1)]: 



































\ X / 


\ e / \ 



(|£;(+)|2 + |£;(-)|2)/2 
im[i;(+) *£;(-)] 



(24) 



The temporal coordinate is the semi-sum of the fluxes at 
each side of the multilayer. The interval remains invari- 
ant 



(e0)2_(ei)2_(e2)2 = if2 



(25) 



This number K (which is the 'radius') is the semi- 
difference of the fluxes at each side of the multilayer; 
i.e., (]£;(+) p - |i?(")p)/4 for both ambient or substrate 
and, therefore, it can take any real value. Without loss 
of generality we can renormalize these variables so as to 
take the value of K equal to 1 and then we are working 
on the unit two-sheeted hypcrboloid of S0(2,l) p^ . 

In summary, given the multilayer matrix Ma.s in 
Eq. (H), the corresponding Lorentz transformation in 
S0(2,l) is I 



HMas) 



2Re(a/3) 
-2Im(a/3) 



2Re(a^*) 
Re(a2 + /32) 

-Im(a2 -I- /32) 



2Im(a/3*) 
Im(a2 _ /32^ 



(26) 

Moreover, using Eq. (^3|) it is direct to obtain the explicit 
expressions in S0(2,l) for each matrix appearing in the 
Iwasawa decomposition (0) . Indeed one flnds 



Ak(<^) 
Aa(0 
An(i') 



1 

cos (j) sin q 
— sin (j) cos (f 

cosh^ — sinh^ 

1 
— sinh ^ cosh f 

V 1 —V 

z/2/2 V l~{v^/2) 



(27) 
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The action of these matrices in S0(2,l) is clear: A^^{(j)) is 



a space rotation of angle (f> in the 



plane, Aa(^) is 



a boost in the direction of the axis with velocity v/c — 
tanh^; and, finally, An(z^) represents a space rotation of 
angle tn [such that tan(TN) — i'/2] followed by a boost 
of angle tn and velocity v/c — tanh(zy/2), both in the 
— plane. In the next Section we shall explore in 
more detail the properties of these three matrices. 



III. GEOMETRICAL INTERPRETATION OF 
THE MULTILAYER ACTION 

A. Multilayer transfer function and hyperbolic 
geometry in the unit disc 

In many instances (e.g., in polarization optics |llo|| ) we 
are interested in the transformation properties of field 
quotients rather than the fields themselves. Therefore, it 
seems natural to consider the complex numbers 



E 



(-) 



e: 



(+) ' 



(28) 



The action of the multilayer given in Eq. (^ can be then 
seen as a function Za = f{zs) that can be appropriately 
called the multilayer transfer function ||l^ . 

From a geometrical viewpoint, this function defines 
a transformation of the complex plane C, mapping the 
point Zs into the point Za, according to 



Za = ^Mas,Zs] 



a + (izs 

and the point of the infinity is given by 



(29) 



$[M,„-a//3] =oo , $[M,„(»] =a7/3 . (30) 

This bilinear transformation defines an action of the 
group SU(1,1) of multilayer transfer matrices on the com- 
plex plane C. The complex plane appears then foliated 
in three regions that remain invariant under the action of 
the group: the unit disc, its boundary and the external 
region. In fact, 



2Re(/3a*z,) 



|a|2 + |/3|2|z,|2 + 2Re(/3a*z,) 



(31) 



and, in consequence, the difference between the numera- 
tor and the denominator is 



{\a\^~m^){\z,\^-l) = \z,[^ 



1 



(32) 



which shows our assertion relative to the invariance of 
those regions. 

Alternatively, the unit disc can be seen as obtained 
from the upper sheet of the unit two-sheet hyperboloid 
defined by Eq. ( |25| ) with K = Ihy means of stereographic 
projection using the south pole (—1,0,0) as projection 
centre. In fact, a simple calculation shows that, in such 



a case, the projection of the point (e", , e^) becomes in 
the complex plane 



1 + e" 



E(+) ' 



(33) 



in full agreement with Eq. 

The boundary of the unit disc corresponds to the pro- 
jection of the infinity point and then, since \za\ — \zs\ = 
1, it can be identified with the action of perfect mirrors 
(i.e., Tas = 0). 

The Iwasawa decomposition has an immediate trans- 
lation in this geometrical framework, and one is led to 
treat separately the action of each one of the matrices 
appearing in this decomposition. To this end, it is worth 
noting that the groups we are considering appear always 
as groups of transformations of some space. The con- 
cept of orbit is especially appropriate for obtaining an 
intuitive meaning of the corresponding action. We recall 
that, given a point P, its orbit is the set of points P' 
obtained from P by the action of all the elements of the 
group. In Fig. 2 we have plotted a typical example of the 
orbits for each one of the subgroups of matrices Ak((^), 
Aa(C), and An(z^). For Ak('?I)) the orbits are the inter- 
section of the hyperboloid with planes ~ constant, 
for Aa(0 with planes = constant, and for A|\i(j/) with 
planes — — constant. 

Through stereographic projection, as indicated in 
Fig. 2, we are working in the unit disc and the corre- 
sponding orbits for the SU(1,1) matrices K, A, and N 
are 



= <i>[K((/)), z] = z exp(-i0) , 

z - i tanh(^/2) 



= <i>[A(e),^] = 



^N{iy),z] 



H-iztanh(^/2) ' 
z + (1 + iz)i^/2 



(34) 



1 + (z - i)i^/2 

As plotted in Fig. 3. a, for matrices K the orbits are cir- 
cumferences centered at the origin. Since for these matri- 
ces Ras = 0, the action of any antireflection system can 
be always pictured as one of these circumferences. For 
the matrices A, they are arcs of circumference centered in 
the real axis and going from the point +i to the point —i 
through z. Finally, for the matrices N the orbits are cir- 
cumferences centered in the imaginary axis and passing 
all of them through the points i, z, and — z*. 

The importance of the Iwasawa decomposition reflects 
also at the geometrical level: no matter how complicated 
a multilayer is, its action can always be viewed in terms of 
these three basic actions with a clear geometric meaning. 
Its explicit application to a real case will be demonstrated 
in the next Section. 



B. Hyperbolic geometry in the upper semi-plane 

The unitary transformation ( p^ has played an impor- 
tant role as intertwining between multilayer and geomet- 
rical optics. One can expect that the structure defined in 
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the unit disc for the former could be translated accord- 
ingly for the latter. To this end, note that if the point 
w e C is defined in terms of z by 



(37) 



z + i 
1 + iz 



(35) 



it is easy to check that the interior of the unit disc is 
mapped onto the upper semi-plane of the complex plane 
w, the boundary maps onto the real axis, while the exte- 
rior of the unit disc becomes the lower semi-plane. 

The relationship between SU(1,1) and SL(2,R) out- 
lined in Section 2 allows us to transport the action of 
SU(1,1) onto C to give an action ^ of SL(2,R) onto C, 
finding in this way an alternative model of the hyperbolic 
geometry and for determining the orbits of the matrices 
of SL(2,R), that are the natural arena of geometrical op- 
tics. 

The corresponding orbits are now 



w — tan(0/2) 



1 + 7«tan(0/2) 
w' - *[^(e),H =«^exp(-C) , 



(36) 



For matrices /C the orbits are circumferences centered in 
the imaginary axis passing through w and —l/w. For 
the matrices A, they are straight lines in the upper semi- 
plane passing through the origin and the point w. Finally, 
for the matrices A/", we have straight lines parallel to the 
real axis passing through the point w. In Fig. 3.b we 
have plotted these orbits in the complex plane w. 



IV. A SIMPLE EXAMPLE AND CONCLUDING 
REMARKS 

It seems pertinent to conclude by showing the power 
of this geometrical approach. In consequence, we shall 
analyze in some detail a practical example: a single glass 
plate of refractive index rii — 1.5 and thickness di — 1.1 
mm embedded in air. The plate is illuminated with a 
monochromatic light of wavelength in vacuo A = 546 nm 
that impinges from both ambient and substrate at an 
angle 6*0 = 45°. 

For this system a standard calculation gives the follow- 
ing reflection and transmission coefHcients: 



roi[l~exp(-i2/3i)] 
1 - rgi exp(-z2/3i) ' 



as 



(1 -r^i)exp(-z/3i) 
exp(— i2/3i) ' 



'01 



where rpi is the Fresnel reflection coefflcient at the inter- 
face 01 (which applies to both p and s polarizations by 
the simple attachment of a subscript p or s) and /3i is the 
plate phase thickness 



2tt 

Pi — —nidi cos6'i. 
A , , 

If we take as initial condition that 



(38) 

the substrate 

Zs — 0.5 exp(i7r/6), then we obtain from Eqs. <m) the 
value Za = —0.6149 -|- 0.0498 i (for s polarization). In 
Fig. 4 we have plotted these points Zg and Za in the unit 
disc. Obviously, from these (experimental) data alone 
we cannot infer at all the possible path for this discrete 
transformation. 

However, the Iwasawa decomposition remedies this se- 
rious drawback: from the geometrical meaning discussed 
before, and once we know the values of (j), ^, and v [that 
are easily computed from Eqs. (^)] we get the intermedi- 
ate values of z' for the ordered application of the matrices 
K((/)), A(^), and N(i/), which, in fact, determines that the 
trajectory from Zs to Za is well defined through the cor- 
responding orbits, as shown in Fig. 4. 

Moreover, and this is the important moral we wish to 
extract from this simple example, if in some experiment 
the values of Zs and Za are measured, one can find, no 
matter how complicated the multilayer is, in a unique 
way, the three arcs of orbits that connect the initial and 
final points in the unit disc. 

We stress that the benefit of this approach lies not in 
any inherent advantage in terms of efficiency in solving 
problems in layered structures. Rather, we expect that 
the formalism presented here could provide a general and 
unifying tool to analyze multilayer performance in an ele- 
gant and concise way that, additionally, is closely related 
to other flelds of physics, such as special relativity and 
geometrical optics. 
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FIG. 1: Wave vectors of the input [Ea and Es ] and output 
[Ei~^ and iJs^'] fields in a multilayer sandwiched between two 
identical semi-infinite ambient and substrate media. 

FIG. 2: Unit hyperboloids defined in Eq. (|2^), which represent 
the space of field states for S0(2,l). In each one of them we 
have plotted a typical orbit for the matrices: a) Ak, b) Aa, 
and c) An . In all the figures we have performed stereographic 
projection from the south pole S of the hyperboloid, to obtain 
the unit disc in the plane = and the corresponding orbits, 
which represent the actions of the SU(1,1) matrices: a) K, b) 
A, and c) N. 



FIG. 3: a) Plot of several orbits in the unit disc of the elements 
of the Iwasawa decomposition K, A, and N for SU(1,1) (from 
left to right, respectively), b) Corresponding orbits in the 
upper complex semiplane for the Iwasawa decomposition /C, 
A, and A/" for SL(2,R). 



9 



FIG. 4: Geometrical representation in the unit disc of the 
action of a single glass plate with the parameters indicated in 
the text. The point Zs is transformed by the plate into the 

point Za- Wc indicate the three orbits given by the Iwasawa 
decomposition and, as a thick line, the trajectory associated 
to the plate action. 



This figure "figure2.gif" is available in "gif" format from: 



http://arXiv.org/ps/physics/0104050vl 
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